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ABSTRACT 


Diieesetnesis developed the ballistics processor of a 
Memenobe Processor airborne tactical system. The multiple 
processor System consisted of three INTEL-8080 microcomputers: 
MiewexeCltive processor, the navigational processor and the 
meets tics processor. the ballistics processor utilized a 
mera! second order Runge-Kutta method of integraticn of 
the equations of motion of unguided air-to-surface weapons. 
imeem ballistics processor computed sufficiently accurate and 
mmely sol@taons to enable the executive processor to extra- 
Mouate an accurate release point for the weapon. The 
fncorithmepermitted complete flexibility in release conditions 
and allowed a complete arsenal of air-to-surface weapons 
mererenitLly canried On the A7-E aircraft. The cost of the 
muersties processor USing current “off the shelf'’ components 
imeplos5 and=the entire tactical system was estimated to ccest 


$3891. 
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ee SUNT RODUCT ION 


fiatary alrborne tactical systems are used to aid the 
mulermand bOombardier/navigator in their specific functions. 
ommeeeattack aircraft, the tactical system is used in the 
following manner: 
feeeO Navigete to a given target 
2. to issue steering commands enroute and/or to guide 
fincwdthe Op 1 LO 
MO COmClc pice mun che Crew the current target location 
imencearron to the aircraft 
weeeto somve ballistics equations which enable the crew 
to release a weapon at the appropriate instant. 
The presently operational systems employ a small 
general purpose computer, such as the IBM 4 PI system. In 
mime enesis, the ballistics processor of a multiple processor 
beeen zation Of a typical tactical system for an attack air- 
Graft is presented. The system's architecture made use of 
three INTEL-8080 microcomputers, each of which was dedicated 
mieecame to itS assigned function. The purpose of the 
Piero le processor organization 1S to exploit the presently 
emerging inexpensive microprocessor technology to solve 
fieeteary Cactical problems. he system's organization 1s 
shown in Figure 1. Because microcomputers presently have 
no hardware multiply-and-divide functions, the programmed 


software mulitiply-and-divide functions caused the microcomputers 
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Meerun On a typical mix of operations approximately six times 
mover than tne presently operational computer systems. By 
meorgdnization of the computations and the use of table 
Meereup techniques, the described system was capable of meet- 
[reece real~time requirements of an airborne weapon delivery 
system. 

The ballistics computer uses the best estimate of the 
fmmeemon, dive angle, altitude, and velocity of the aircraft 
Mts COMPutation. It is ajso provided with the weapon 
miemertication (for the particular drag characteristics) 
and the wind velocities. From this information and the 
target's altitude (above Mean Sea Level), the computer 
determines the down range travel and time cof fall of the 
mm@ouaded weapon (1.e., the impact coordinates in the target 
plane). 

The accuracy with which a ballistic weapcen can be 
delivered against a target depends greatly upon the accuracy 
of the sensor-supplied release parameters as weil as the 
meeuracy with which the ballistics equations are solved to 
Meedict the weapon's down range travel. 

Most aircraft unguided air-to-ground weapons can be 
Seeeribea as ballistic projectiles. The only forces acting 
on them after release from the aircraft are gravity and 
aerodynamic drag. Bullets, streamlinec bombs, dragued 
(retarded) bombs, cluster munitions and unguided rockets 


PPeeeremovormyurmout) are all ballistic projectiles. Guided 


eZ 





Meapons and weapons developing lift are not ballistic pro- 
Meewiles nor are they considered in this development. 

This thesis first describes the over-all tactical system 
mumorder to give the reader a perspective of how the ballistics 
processor functions within the entire weapon system. Included 
metnis background material is a discussion of the naviga- 
mar computer, the ballistics computer, the executive com- 
Smeer and the intcrcomputer communications network. 

mie Dasic Calculations involved in solving the ballistics 
Preeorem 1S described next. The assumptions used in the pro- 
mime tOrmulation are stated; the derivation of the particular 
differential equations is presented; and the application of 
the Runge-Kutta integration method is shown to be suitable 
for current operational airborne weapon systems. The four 
general classes of weapons and how the algorithm handles eacn 
/@ficeerdtion is defined. The unique method of compressing the 
drag data for all twenty-eight weapons is shown to be a simple 
curve fitting of the drag coefficients with respect to scme 
Standard. 

The simulation was initially programmed using FORTRAN 
for ease of testing and debugging. The important phases in 
programming were in the DERIV and RUNGE subroutines where a 
Slose Study of the optimum calculation methods was conducted. 
The program was tested for accuracy of results by comparison ; 
with the published range tables (NAVAIR 01-1C-1T-1). e 

Upon completion of the FORTRAN implementation, the re- 


Subbale Program was translated into PL/M, a high level 


IES 


A 





Janguage compatible with the selected microcomputer (INTELLEC 
Be). The selected test cases were verified correct in the 


PL/M version of the program. 
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II, BACKGROUND 


meorder tO give the reader a complete understanding of 
memetne ballistics processor functions, this section is 
fevered to describing the overall tactical system in which 
Maemballistics processor is an important part. The first 
section describes the function of the navigational computer, 
the sensors which provide data to the computer and the 
meeorrthm used to compute present position. 

ites bpallaistics computer will be completely defined in 
miemremaining sections of this thesis and will only be 
Metetly reviewed here. Using the best estimate of position 
and velocity at a given instant, together with information 
meeeene weapon's drag characteristics, the ballistics 
computer solves a system of four differential equations and 
@erculates the impact coordinates of the weapon in the target 
plane. 

Witemercctiiive Computer controls the dispiay and generates 
commands to the autopilot and weapon release mechanism. This 
memoucter iS interrupted asynchronously as fresh data are 
veloped in the navigational and the ballistics computers. 
iiiewexecutive computer extrapolates from the data the appro- 
priate release instant of the weapon. 

Wicmiieerecemputer communication net is described next. 
The star-like communication net provided a simple and rapid 


Means Of ianformation transfer. 





Pee tHe NAVIGATIONAL COMPUTER 

This subsystem 1S present in all tactical systems. In 
operational systems, a subprogram (module) is periodically 
executed to update the present position by the change in 
Meertion in the last time increment. 

hmenesscompucer, the sensor instruments provide digital 
meee to the cemputer at specified time intervals, controlled 
Meeeoeclock. tIhere are four sets of sensor instruments which 
provide the data: 

1. The Inertial Navigation System (INS) 

Pee ihe Doppler Radar (DR) 

3. Air Data (AD) 

meee Radar Set (RS). 
The INS system provides attitude (roll and elevation) and 
fmoerty increments in the x, y, z directions. The DR is a 
Mere tty sensor that uses the doppler principle for the 
continuous measurement of aircraft ground track speed and 
Meert angle. The AD computes outputs of modified corrected 
Static pressure, pressure altitude and mach number from 
meiete anid pi1tot tube pressure inputs. The RS provides 
mareect azimuth, elevation, and range signals to the naviga- 
meenal Computer. A magnetic compass measures the aircraft 
heading in relation to magnetic north. In order to navigate 
Meeurately ian this mode the wind speed and direction are 
g@etermined by using the “sensed" velocities of the aircraft. 
Altimeter measurements are also sampled in this system to 


update the present position altitude. 
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iene are four basic modes in which this system operates. 
Mme crew has a choice of using both the INS {velocities} and 
emer eputs tO minimize the Schuler oscillations. This is the 
most reliable and accurate mode. If the crew suspects that 
eueemer the INS or the DR system is malfunctioning, they may 
use only the data from one of the sources. The least 
reliable and inaccurate mode is the last resort, when both 
the DR and the INS are malfunctioning. The AD mode uses the 
meaeeriially estimated wind vector and true airspeed of the air- 
meri tO determine the position in latitude and longitude. 

miedibeoperatimeg modes, the same basic technique of smooth- 
mieminput data is used. The least squares curve fitting 
meemnague 1S used to fit a quadratic polynomial to the most 
recent n (n =< 16} data values. Use of the Legendre 
Seenogonal polynomial functions is made in order to make the 
numerical calculations as accurate and rapid as possible. 
mmoouhing the short time increment data values will generate 
Meese total error in the results. 

Msane the smoothed data values, the appropriate calcu- 
Mites are made to estimate the distance increments in the 
Memeerime Interval. As soon as the new value is generated, 
eieedrsplay of the navigational computer is refreshed, the 
Secucive Computer is interrupted and this most recent 
eeeeamate iS transferred to the executive computer. If the 
navigational computer is functional, the crew has informaticn 
Seeeene Present position, independent of the executive 


computer. 





The least squares curve fitting technique enables the 
System to diagnose sensor malfunction quite easily. If the 
sensor generates fluctuations in the input data above a 
certain threshold, then the sensor can be diagnosed as mal- 
memersoning. Similarly, if two sensors disagree above an 
allowable threshold, then an automatic data analysis will 
Mmeeest which senscr is failing. If the crew concurs in 
this diagnosis, the operating mode of the navigational 


meeeem can be changed to exclude the faulty sensor. 


Pine BALLISTICS COMPUTER 

maemo tbhisties computer 1s provided with the present 
einates Of the positaon, dive angle, altitude, and 
memoecity of the aircraft. It is also provided with the type 
of weapon and wind conditions. From this information and 
Mm@emtarget's altitude, the computer determines the down- 
maee travel of the weapon, if released at this instant. 
The assumption of no guidance after release from the air- 
craft is still applicable. However, the weapon may have a 


rocket assist, or a retarded mode of fall. 


Meise EXECUTIVE COMPUTER 

The executive computer controls the main display and 
Pemerates Steering commands to the autopilot and firing 
peses tO the weapons release mechanisms. One of its 
Peeetecal Lunctions is to extrapolate from the data generated 


by the ballistics computer the weapon release point. 
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Ec cimoMmlencloc. 1S used to record the time at which 
the asynchronous data were generated by the ballistics 
empmter, Ine position of the impact point in the target 
mumemencan be described by a quadratic function of time. 
Miererore, the real time release point (adjusted for time 
mommays) Can be extrapolated very accurately from these data 


Parenout additional extensive calculations. 


Petits INTERCOMPUTER COMMUNICATIONS NET 

In Order to keep the design of the intercomputer 
mommunications net simple and inexpensive, the starlike 
network was chosen. The executive computer is the center of 
the star and one input port is multiplexed to the two output 
ports of the navigational and ballistics computer as shown 
Mierigure 2. 

The intercommunication takes place in the following 
Mapner: An initializing phase of the program interrogates 
mre keyboard to obtain data about the initial position of 
Meemeatrerait. Ihe latitude and longitude are entered in 
degrees and minutes to the nearest 1/10 minute; the 
altitude is initialized at the height above mean sea level. 
These data are made available to the navigational computer. 
Also the type of weapon and the position of the target (if 
Known in advance) is made available to the ballistics 
Computer. Similarly, the mode of navigation is determined 
(if possible prior to flight) and relayed to the navigational 


Sompucer. 
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The execution phase of the computer program makes use of 
miewreal time clock in the executive computer which is 
feym@enronously interruping the navigational computer in 
emegenre that sensor information be transferred via input port 2 
to the navigational computer. When the navigational computer 
feeecompleted its calculations, it interrupts the executive 
mememcer and transfers the calculated current position 
Begmentidally via input port 1 to the main computer. The 
ameunt Of information sent is 9 bytes, and the transfer is 
@eeomplaished in about 250 micro-seconds. 

Similarly, when the ballistics computer has completed 
memmeatculations, it interrupts the executive computer and 
meemsrers its information sequentially via input port 1 to 
Mmemexecutive computer. Only 6 bytes of information are 
Semen approximately 170 micro-seconds. Between interrup- 
tions, the executive computer updates the display using 
output port 2. Calculations are continuously made to extra- 
polate the release instant of the weapon, and the release 


mechanism is activated when appropriate via output port 3. 





Peo nkololles PRODGEM 


aS A ee et 


The equations of motion governing the trajectory of a 
maeenstic projectile are a simple-appearing set of second- 
Order differential equations whenever the drag is neglected. 
However, when a reasonably accurate model of the aerodynamic 
freeecaused by the projectile's motion through the air mass 
gS included in these equations they are rendered nonlinear, 
and no general] solution in closed form has been found that 
Mmemsacistactory for solution by an airborne fire-control 


computer. 


fee BASIC ASSUMPTIONS 

The mathematical model was chosen to represent reality 
as accurately as possible. The armed forces publish range 
tables for various weapons and therefore the same model was 
chosen as the one used to construct these tables (NAVAIR 01- 
me-jT-1). 

The equations of motion were developed assuming the 
meogectilé 1S a point mass acted on only by the force of 
Meavity and the retardation forces due to air resistance. 
The trajectory was restricted to a piane by ignoring the 
Beesstrail effects of the wind. These crosstrail effects 
ana the component of wind which affects drag are considered 


mieene executive computer. 


(ab 





The assumptions are summarized below: 

feline Barth 1s flat and nonrotating. 

fee ihe gravitational attraction 1s constant. 

fee ihe projectile 1S a point mass. 

eine projectile is not powered (except rockets) and 


mas a, constant mas®?. 


Pee DERIVATION OF THE EQUATIONS 

The weapon's fall is assumed to be governed by a drag 
meee which is proportional to the square of the velocity. 
Newton's equation of motion becomes 


Hien ty dt = = 8 Mm -|Ivel| Cw vt/ || ve || 


where the only two forces acting on the weapon are gravity 


and aerodynamic drag. 


tie total Velocity 
Mmeceiiass OL the projectile 
e - force due to the Earth's gravitation 
five; drae coefficient which depends on the 
characteristics of the weapon and may vary 
with time (in case the weapon has two drag 
elie es)\ and it Varies Slightly with airspeed. 
meerder to solve this differential equation, initial 
monaytions, that is the position and velocity, must be 
mgecttied, the most convenient coordinate system is height 
me wanpove sea ievel and "x" downrange travel, which is zero 


at the time of release. 





In the x and y component form, Newton's equation becomes 


dVx/dt = -{{ve || /m cove: (Cw 

dVy/dt = - g - |\vell ae ve OW. 

dx/dt = Vx 

Dy/dt = Vy 

x (0) = 0 downrange travel 

CO) = present altitude 

Vs(0) = present horizontal velocity 
Vy(0) = present vertical velocity 


meee RUNGE-KUTTA APPLICATION 

Maeechoice of the “best' Runge-Kutta (R-K) formula for 
Meer ticular problem was not entirely obvious. The lower 
Order formulae are simple, but require more integration 
Seems than the more accurate higher order formulae. On the 
Gener hand, the complexity of the higher order formulae 
Teeeereguire so much computer time that this becomes the 
Overriding consideration. Generally, it was found to be 
Safer to use the lower order formulae when discontinuities 
in the coefficients of the differential equations can be 
expected. Such discontinuities frequently occur in bomb 
mempectories (e€.¢., drogue deployment or thrust termination). 
The second-order R-K formula also afforded the greatest 
Meedbility in the distribution and size of the integration 
meee. these considerations, along with their simplicity, 
led to the choice of the second-order R-K method as the most 


Spee epriate for the ballistics algorithm. 
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fei) = Wx Fl = -1/m* |[vt|| * vx * Cw 
U2 = y te e= VY 
U3 = Vy F4 = -1/m *|/vt|| * vy - cw - g 


The differential equations may be expressed as 
Ai dee bails 2503, 04 ) shee a ee | 


Therefore, the Runge-Kutta second order method used in the 


program becomes 


mli = At * Fi(U1,U2,U3,U4) 1 e384 
m2i = At * Fi(U1 + cAtm11,U2 + cAtm12,U3 + cAtm13, 


U4 + cAtmi4) Oe = 942,354 


+ 


a. Ue ly omc ecre-wrije Mliet m21) a = 1,2,3,4 


The procedure computes new values for the variables 
meer. the remaining time of fall, DIV, (as computed for a 
vacuum release), is less than the value of the next integra- 
tion time increment. The procedure then uses the value of 
Meetor the final integration step size, and in the calcula- 


tion of the down range travel of the weapon. 


Meee CLASSES OF WEAPONS 
Meeoencamlined Bombs and Bullets 
The process described above represents streamlined 
bombs quite well. After being fired, a bullet is really just 
a small, streamlined bomb and was handled as such. The muzzle 
fooecity (5500 Lcet per second) of a particular gun-bullet 


combination was added to the aircraft velocity. 
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2. Drogued Bombs 
These@pombs are more compiex in their trajectory 
modeling. Their common characteristic is that they are all 
released in a relatively low-drag configuration. At some 
time after release, they deploy vanes or a parachute which 
mmeatly increases the totait drag. To compensate for this 
mmeamee Of drag characteristics, the initial integration step 
Size was chosen to coincide with the deployment of the high 
Mmeeomdeyices. Thereafter, there was no difficulty in calcu- 
Beene the trajectory of the drogued weapons. 
fee Ciuster Bombs 
These weapons are released as large, low-drag 
Mmemmeainers. At a predetermined point in their trajectory, 
Which is time dependent, they use some method to dispense 
iamber weapons. These smaller weapons are of higher drag 
Mm@am the container. The ballistics computer will predict 
miemimpact point of the center of the pattern rather than 
individual positions of the bomblets. The method of applying 
meemoieOrithm 1S to calculate the trajectory of the container 
memeene POInt where it dispenses the small weapons; switch 
mea ditierent drag curve which describes the motion of the 
Mmmeeern center; and integrate down to the target altitude. 
As in the retarded bombs, the initial integration step was 
chosen to coincide with the act of dispensing the small 
weapons. 
poe lcneuided Rockets 
POGkets present Several unique complexities. They 


Change weight during flight; they have thrust as well as 
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feete, and they slew around aiter firing, because their 
mmenes May not line up with the aircraft direction of 
Mmenpent. ihe algorithm treats unguided rockets like a 
@eister bomb with the first stage having an average thrust 
@oewell as an average thrust time. This provides sufficient 
accuracy and saves computer storage that would be taken 


me@ieine the thrust profile of each rocket. 


E. COMPRESSION OF DRAG DATA 

Mieebpechavior of a particular weapon's drag coefficient, 
Cd, as a function of Mach number must be known to calculate 
impact ranges for that weapon. Normally, this information 
is given in tabular form with some weapons having more than 
meeepoints in the table. Most aircraft have the capability 
memcarry many different types of weapons. It is desirable 
memestore all the Cd information of an aircraft's weapon 
mpertoire in the airborne computer to avoid loading the Cd 
data for each different combination of weapons. This 
Mmeaguares Storing drag coefficient data in a more efficient 
form than a table. 

The following factors influence the choice of a scheme 
tO approximate the drag coefficients for a given set of 
weapons: 

I. Computer storage 

2. Computer computation time 

5. Delivery envelope of each weapon 


4. Allowable downrange and time-of-fall errors. 


27 





ihemrould be ideal to use one general, fairly powerful 
expression to fit all the drag tables for the weapon 
Memfertoire of a given aircraft. This is quite difficult to 
accomplish because weapons like guns and rockets require 
fmmveiues for Mach numbers much higher than free-fall bombs. 
Pemeeeconsequence of this generally acceptable tabling pro- 
Soaure, extra storage is spent on weapons that do not require 
Mm@emextra capability. 

mesroniticant Saving can be realized when it is 
recognized that several free-fall weapons have drag 
coefficient curves that ditfer only by a multiplicative 
m@enmor. the Mk 80 low-drag series are one such group of 
Weapons with the Mk 84 drag curve conventially taken as the 
reference. Similarly, the Mk 106, CBU bomblet, and Sadeye 
pomelet are another group with the Garve drag curve as the 
mmeemence. The specific region of interest of many more 
drag curves can be approximated by multiplying these two 
Mmemmerence drag curves by a factor and then translating them 
meome the Cd and M axes until they match the original drag 
SUnVves . 

To further aid in the explanation of approximating 
weapon drag curves by this method, let Cd(M) and M be the 
drag coefficient and Mach number of a weapon drag curve. 
iemeeciM) be the drag coefficient of a reference drag curve 
eaaetet A, B and DM be the multiplicative factor, drag 
Seetticient translation, and Mach number translation,respec- 


tively, that when applied to the reference drag curve wil} 
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elt 


“proximate the particular weapon drag curve of interest. 


That is: 
Cd(M) # A # Cd(M + DM) + B 


A hypothetical weapon drag curve, Cd, and reference drag 
meen CG, are shown in Figure 3. For simplicity of discus- 
Sion, the Mach number intervals (M2 -M1) and (M4 - M3) are 
eee. Suppose the Cd curve is multiplied by an A so that 
Meemeceand the A * Cd curves are nearly identical. For the 
mmemoann ~ Ca curves to be a congruent, the A * Cd curve must 
be translated along the Mach axis an amount equal to DM, and 
meee the Drag Coefficient axis an amount equal to B. 

For many weapons, the above stated equation will not give 
an accurate fit over a wide interva] of Mach numbers; instead, 
Eeeeeanc DM snould be chosen so that the Cd fit is best for 
Siem ach number interval that corresponds to the delivery en- 
femepe fOr a particular aircraft. 

Considering the weapons currently implemented on a typical 
attack aircraft, (A7-E), both the Mk 84 and the Garve reference 
drag curves are necessary to handle the free-fall weapons. 
Instead of extending these curves to handle guns and rockets 
that have a much larger Mach number range, a rocket-gun 


Memerence drag curve is used. 


PeecURVE FITTING THE DRAG COEFFICIENTS 
In the previous section, the reference curves needed and 


how they were used to approximate weapon drag curves for use 
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an the algorithm was presented. In this section the reference 
Mmereecurves versus Mach number is represented. The form of 
emeyrunction used will minimize computer storage, logic and 
computation time for the allowable error in downrange and 

time of fall for the different weapons. It was decided that 
Gach of the three reference curves be divided into three 
meootens With Cad expressed as a second-order polynomial in 


Mach number for each region. 


Cd(M) = a0 + al * M+ a2 * M? 


Mimewcocificients, aO, al, and aZ depend on the region the 
Mach number M is in. The drag coefficient for a particular 


weapon is expressed as: 
Cd(M) * A * (a0 + al * (M+ DM) + a2 * (M+ DM)*) +B 


Meeea tUrther step to save computer storage and reduce the 
computation time for the algorithm, another compression of 
meee data was undertaken. The computation C * (PI/8 * Cd), 
MeG, is made in the Runge-Kutta portion of the algorithm 
where C is the bomb factor, aie where d is the weapon 
fiemeter in feet and W is the mass in slugs. One constant 
per weapon in computer storage and two multiplications per 


integration step is saved by expressing CKDG as 


CKDG 


ie echoG a elNG 


where CE = A * (d/d)* * W/W 





CKDG = d“/W * PI/8 * Cd 
DKG = d7/W * PI/8 * B 


ieee above equations, d and W are the assumed diameter 
memmdss for the reference curve Cd. The final form of CKDG 
in the algorithm is 


CKDG(M) = CF * (bO + bl * (M+ DM) + b2 * (M + DM)“) + DKG 


mete | gives a list of bO, bi, and b2 for each region of the 
three reference curves and the values of CF, DM, and DKG, are 
tabulated in Table 2. 

Mmiimust appears in the algorithm in the total drag 


function 


Bee -sin/V¥ - RHO“ CKDG * V 


where TH tamist,mass 


RHO = air density 


V = weapon velocity 


mmomckKDG was previously defined in the above equations. 
Heeure 4 shows thrust versus time for a typical rocket 
approximated by a constant, thrust = total impulse/Tl, 

where Tl is decreased from actual motor burn time in order 
to make thrust equal to the average thrust of region 2. The 
mass of the rocket, W, obviously decreases during rocket 


motor burn time. However, it is adequate for W to be 
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approximated by a constant, by varying W until downrange 
error is minimized. The second stage of the rocket, for time 
meer weapon release greater than Tl, is free fall and TH 


will be zero. 
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PV. FORTRAN SIMULATION 


feeebOMB DROP SIMULATION 

Mero initial ballistic trajectory algorithm for digital 
airborne fire control computers used in the simulation of 
Traj, Deriv, and the Runge routines, was developed by the 
Boeing Company, Seattle, Washington. The research was 
Sempketed under a Naval Weapons Center contract in September 
970. 

AS previously mentioned, the Second-Order Runge-Kutta 
integration scheme was chosen because of its flexibility 
and simplicity. The results, as evidenced in Tables 3, 4 
and 5, compared favorably to the published data (NAVAIR 
ieee -lT-1), with mean errors of 0.10% for the down range 
travel. The low-drag bombs shoved relatively little 
Sensitivity to the step size, so long as the maximum 
mitegration step, (DMAX), was less than or equal to 5 
seconds and the normal Initial Integration Step, (DTI), 
was less than or equal to 3 seconds. However, the dual 
Stage weapons (IDNO>17) did show a large sensitivity to 
the step size used during the integration scheme. The 
Peeuracy desired, required the DTI value to be consistently 
less than or equai to 2 seconds with the average value of 
1.05 seconds, for the final eleven weapons. 

The ultimate verification of the step size for the 


mec -hutta ({Meriv also) method was conducted by varying 





meres. oample Errors in Impact Range Caiculations for the 
General-Purpose Bomb Mk 82/Snakeye/Unretarded. 









Range fi fImpact Range 


Milgorithm error, 





Reblease Conditions mmia.c.t 


MeLocity, Aogeke .i ballistic 
knots iE ie deg | Tables _ 
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400 2 5 B10 0 oO) Si 0 
5000 = 8, 9153 steal 
LOO = 2110 FOS 36 rol 
warere -45 on 0 = | 
450 ZO 0 Gi 3 =f 
5000 =10 reo2> a 
12000 =o0 11340 +-] 
2000 -45 7298 te 
500 Zo 0 0 9946 9940) 1 = 6 
5000 a) 10843 1S 3.9 -4 
12000 = 5), P20 ov 7067 0 
12000 -45 8406 8407 rl 


fee 4. Sample Errors in Impact Range Calculations for the 
Special-Purpose Bomb Mk 40/Mcd O/Conical Tail. 
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mooie Do. sample Errors in Impact Range Calculations for the 
Practise-Purpose Bomb Mk 76/ with Lug. 


Release Conditions impact Range, Otis impact Range 


memmoecity, Altitude,|[Angle,}| Ballistic fAlgorithm| error 
knots ae deg Bab bers 


400 2500 0 718) 318 TON 7 = Ih 
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12000 ~45 ee W235 ae 
450 2500 0 8463 84.99 -4 
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12000 -45 7669 7674 17S 
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the time interval from one to five seconds. The TRAJ 
Hicegration scheme was simultaneously tested using several 
Metease CONditions to demonstrate the model's flexibility in 
Simulating bomb loft, level flight release and dive bombing 
conditions. The emphasis, at this point, was on the number 
Mme cional calls to RUNGE for the projectile to reach the 
Mmereet altitude. The results, by altering At, showed changes 
mone from 0.003 to 0.44% for the time of fall, and 0.00% 

to 0.46% for the downrange travel, with a substantial decrease 


fm_meneenumber of cails to RUNGE also evidenced (Table 6). 


feo veRALL SIMULATION 

The PRINT routine was a means of simulating the various 
Seeweit displays in the tactical attack aircraft. The output 
consisted of the release altitude, dive angle, and airspeed, 
as well as the resulting downrange travel and the time of 
fall of the weapon. 

The INPUT routine simulated the collecting and passing 
Seeeenecs Sensor data to the various processors. The sensors 
merece d 2 the attack aircraft are the Inertial Navigation 
System, Doppler Radar, Air Data and the Radar Set. Specifi- 
Cally, the velocity, dive angle, altitude and the identifi- 
cation number of the weapon (usually provided via the armament 
panel) are required for the execution of the program. 

mirc SEIDAT routine completed the interface between the 
Seeeecutive computer and the ballistics computer. Its purpose 


was that of assigning the constants their specific nonvarying 


4] 





ORE TRU: Ea SA ae NS ge A i. 


FSeTiZtoO = SSeTTl cio Osco =) oe c 1a 








Loc ie Ones 










Z00°0 13 O00L 
Ue), (7) : (¢) | (2) : (77) S34 OSS 
ad RE GG ee REN SS es Gy T= 


SHEN ZL0°- 67°C ~06- 














































=e 09641770 =" 8 6964700 = 6.0967: 710 8on TOG 710 = ole! a3 OOS 
Ca.) (7) (¢) IL), : Cou) : S24 00 
TO Ceo Oe Cigi200) 0 | NS acl) oma me scales C10 amo mnan BO) A 
Z7O'- ITZTeuzto’- L110 = “Gren? O° — Soler eVEy 00 (0) 
(¢) (6) (8T) 3 S40 7 
Zit SOs i 0 





O 











LGTS 1 











i,¢0 Se See or CS 1G yO 35 OOQST 
ik) | (9T) Ce SIM OOE 
£608 = 00 


BOO TO) GO 18c2) 
Ce) 





Sve) 
ec” = cena 






MEP — PG RER 










































(jee “Ser Oclvec>— soreodsel = PAG ON 23 OOOE 
(cc) (9) | Che) G8°0Z S24 069 
Mei G61 = 18 0c (Lez %70°- 78°07I | O01 
POR” = OAs = Meee Pet, lek yal +3 005 
(¢) : | (oT) 68°€T S44 0S9 
oO C1 | 4GG = 7e7e 0 1 j4vt'- LEE] ZL0°- Bere o OT 
| Bo ‘ f : meas ole 
= S186 |AG02= 20 G8G4C0°= EB eeeG | SO8— Bene wigs Gees 9786 1s Oe 
(7) ! (¢) (9) | (6) Comp) i 6G aa S34 OSE 
= GIN ESO eT era Be LOO Ee om) 6B ORS non 
Hele= 116699200 - 2 6699X%60°— 8 °8699|7%z20-— ae 7TO0°- 66694 T OH9 23 OO00T 3 
| (€) : (€) (7) (9) (ZT) iL it Ae 'S2 moe 
Sealine Ouaeli ele 60) 1 yal ClO eee Oia aeiaia UOMO) (Cage dete : oOT 
——— : . ; ‘2 pore 5 , , — ; i eee 





OG eS ne ee = nya O°'€ = la | | xSTTAVE SNoraianoo | 
iS AS ASVATIY | 


(€8 UN ‘Teazequy; ewT2) sosey 2S9y ueAIOg ‘9 ZIEVI 





42 





[ge COL RE TE 
Yd Qnd-Yd 3seL aZuey uMog 


(UOTIBASARVI asuny 
Cae S ya Co gf O ON ) 


Ld dnd 


Id GNd-1d 38°] TTes JO ouTy AMLNA ATANVS 









L=bl-st1 0vtt AUN 


SE FT aT PS IT YS SES. (RY FONE OYE ES 
OP ELLE te 


LZO'- T°€6THZZO"- ZEGTEL ZITO’ - Y'ESTEZIO'- YESTEt V6TE | | J5 OOOL 


















On ae 61 E) 






5) ) E50) | 7 (OT) 19°6 S34 00% 
ply = eee S 61201 7O— 09 a 2Ol = O96 09°6 Z00°O 19°6 | Reels 


















33 000 
S34 0S4 
507- 


AU Om aos) 7 20 Oe = imas 9 
(2) Ce) 


%00°O *689% 
OS) 


%0°O 1°689% 689% 
(6) | SI's 
OO Olan uaa 






ae ISG £1 ° 8] 










xSHTAVI 
D1ES 1 ee ol 





SNOT Taos 
A Savedeel a a 


Oe eG 


(penuyf _uos) seseDp YSo9] ueayrzAog °9 AIAVIL 


43 





fees and reinitializing the variables to zero at the start 
of each new weapon's integration. The routine also performed 
various calculations including the determination of the 
mmieetty in the x and y directions, the conversion of raw 
velocity to feet per second and the calculation of the release 
angle in radians. 

The weapon dependent constants required for the algorithm 
were stored in DECODE. Also, DECODE contains the necessary 
Gioeeae tO Load the correct constants for the particular 
weapon the crew has selected. After the first pass, there 
memo need to go through DECODE before entering the algorithm 
memeecrajectory calculation. The actual scheme used for 
SeeoreE in reality will depend on the particular aircraft ana 
ats fire control requirements and how the algorithm is inter- 


imeteed in the airborne computer. 
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V5 PRY MOIMPLEMENTATION 


foeimplement the model on the INTEL 8080 microcomputer, 
the floating point arithmetic package was used to maintain a 
four nexadecimal digit accuracy during calculations. The 
floating point variable therefore occupies three eight-bit 
words: two words for the variable mantissa and one word for 
the mantissa sign bit and the variable exponent. The man- 
tissa of the variable is viewed as a normalized binary number 
With the left most digit always one. 

In order to manipulate these variables, an arithmetic 
Smieroutine library was supplied by Professor Kodres. At 
Mmesent, there exist routines that add, subtract, multiply, 
ftde and take the square roots of these floating point 
variables. 

The actual PL/M version of the ballistics program was 
translated directly from the working FORTRAN version. All 
the associated flow charts, process graphs, program listings, 
and the simulated cutput can be found in the appendix. 

The PL/M routines have been compiled and tested. The 
results compare very favorably with the FORTRAN results. 

The PL/M verification of the test parameters is included in 


Maple 7. 
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Meee 7. PL/M lest Verification (time interval, MK 83) 
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VI. CONCLUSIONS 


The objective of this thesis was to demonstrate a micro- 
computer based digital system which will perform the same 
functions as the attack portion of a presently implemented 
mmeporne tactical system. The main purpose of the project 
was to verify the feasibility of implementing the most com- 
Mlexesarithmetic task of the airborne tactical computer by 
use of an INTEL 8080 microcomputer at an estimated cost of 
$3891 for the entire tactical system. 

Because of the rapidly decreasing hardware costs of the 
microcomputer and significant increases in processing speeds 
Which the microcomputer industry 1s presently realizing, 
tec NeW, innovative systems, which are in all ways equivalent 
to the presently available avionics systems, can be developed 
@t)a hardware cost which is nearly 10 times smaller. The 
sottware development costs are also likely to be smaller 
Because the system can be subdivided into smaller, independ- 
ently operating subsystems, each of which will perform a well 
@eained task. 

The lack of hardware multiply-and-divide operations in 
the microcomputer is compensated by careful process analysis 
memeneauce the number of multiply-and-divide operations. This 
allows three “slow (in the computer sense) computers to 
meee tie same problem that one “fast'' computer is presently 
mmeyine, at a creat hardware cost reduction. The necessity 


for careful analysis of the problem may even disappear soon, 
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because the microcomputer manufacturers are likely to 
incorporate the hardware multiply-and-divide operations. 
This thesis demonstrates the feasibility of using the micro- 
Selputer as a systems’ building block for airborne tactical 


systems. 





Per eNDOT PROGRAM FLOW CHARTS AND PROCESS GRAPHS 
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eymools Used in Process Graphs _ 
pmeyelens output 
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MAIN PROGRAM 


START 


“ll 


SET =6= 
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I DNO=NUMID(IIDNO) 


SET = 0.0 
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MAIN PROGRAM (PART 2) 
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MAIN PROGRAM (PART 3) ey 


SW ETCH 


V/ 








ee) 





mPUL SUBROUTINE 
JIN Seg | 


WAS E[GKG NL Nie 


V/ 


Me EMO) Bb 


V 


DIVE ANGLE 





WEAPON ID. 
NUMBER 


ews 


a. 


54 





PRINT SUBROUTINE 





ENTRY 


V 






GLOBAL 






VARIABLES 





\/ 


(0) 
Vv 


V/ 


RETURN 


DS 






HEADING 







INFORMATION 





SETDAT SUBROUTINE 
ENTRY 


/ 


GLOBAL 









VARIABLES 
G= 32.174 
RAD= .0174533 (0) 
V 
CFORM1= 0.0 
CFORM2= 0.0 
DM1= 0.0 
DM2= 0.0 
DKG1= 0.0 
DKG2= 0.0 
VMUZ= 0.0 
VE= 0.0 
SL= 0.0 
DMAX= 5.0 
[***ITYPE= -] 
FN= 0.0 
IBOTH= 1 
**T= 0.0 
RETURN 


SeohT VALUES FOR PL/M ARE 1, 2, 3 


meet iS @M AND A IS mS IN PL/M. 


B= > iTYPE HAS A VALUE 3 FOR -1 IN PL/M. 
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U= .VKTS (IV) *1.6878 

DEL= ATAN(VE,U) 

V=SORT (U*U+VE* VE) 

VXA= (V+VMUZ) *COS (THETA-DEL) 
VYZ= (V+VMUZ) *SIN(THETA-DEL) 





















DECODE SUBROUTINE a 
ENTRY 


GLOBAL 


VARIABLES 





IREF= 4 
DNGl—= Soe 
bite ao 








TREF= 4 
DKG1= 6.29E-03 
DTI= 3. 


(:)- “{IREF= 4 


DKGl=—4501E- 03 















TREF= 2 
| pMAX= 3. 
2 











CFORM1=.3.92E-03 
DKG1= yone Os 
[DTI= 2. 
(s)}— IREF= 2 
DMAX= 3. 
CFORMI= 3.91E-03 
DKGl= 6.36E-03 


+) 
4 e 


DTI=_ 


DMAX= 2. 
DKGI=" 02021266 
DTI= 1. 


=a a :) 
<a CoO G 
OYI= 3. : 




















DECODE SUBROUTINE 


CEA ee 2 ) 





IREF= 4 
DMAX = 3. 


CFORMI= 2.064 


j= 
DLs. 




























q 
aes 


CFORM1= 1.4932 
DIt= 2. 











kREF= 
CBORM] = 1.3431 
DII= 1. 








TREP= 1 
CFORMI= 1.21 
DTI= 3. 









PRER— 1) 
Cron) 3.12 
BKei—s- te 22E-03 
DTI= 3. 




















TREF= 1 
DMAX= 3. 
CFORM1= 3.4972 
Bi 2 : 






EREF= 1 
CFORMI= 1.605 
Dik= =3- 








— IREF= 4 — >) 
DMAX= 3. 
DKG1= 7.33F-03 
DTI= 1. 
+ vqntmmrmnnmmanienetnateny wancemnamaraniaadh, 











meeODs SUBROUTINE (PART 3) 












fryer 
IBOTH= 2 

DKGI= /.33E-03 
CrhO RZ — sik. 69k —- Ur 


foM2=— 0.38 
DKG2Z= 0.17166 
DS= 0.6617 
siu= OO 000267 


DII= 2. 


TREF= | 
ITYPE= l 
IBOTH= 2 
GEORMI= 2.0/54 
wmerORM2= 0.2217 
DS= 4.267 

DET yl. 5 










IREF= 1 
Liye E= 1 
IBOTH= 2 
ChORMI= 2.2973 
CM1= 0.32 
BRGI= 76-17-03 
CFORM2= 1.11E-02| 
DM2=0.41 
DKG2=0.16885 
DS= 4.06 
DTI= 2.0 






















REE = 1 
IMB G2)9 owe 
IBOTH= 2 
CFORMI= 2.2404 
CFORM2Z—. 0.1175 
DS= 4.0 
Dira— 1262 


rss ) 
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DECODE SUBROUTINE 


(PART 4) 








Pea e= 1 
ITYPE= 1 

IBOTH= 2 

DiGloy9. 7/7E=-03 
CFORMZ— 2. 31E=-02 











DM2= 0.38 
DKG2= 0.23287 
DS= 0.679 
SsL= -0.000303 
Dit— wi 622 










IREF= 3 
DMAX= 1.5 
CFORM1= 2.9964 
DKG1l= -0.014992 
VMUZ= 3300. 
DTI= 0.5 













IREF= 3 
ITYPE= 2 
GFORMI—— 0. 32 
CFORM2= 1.0 
FN= 1746. 
DS —eia4225 
Dri= 1). 
















IREF= 4 
ITYPE= 0 
DS= 0.968 
DKG2= 1.48 
Dri=- O31 











IREF= 4 
ILS) ae 
DS= 0.89 
DKG2= 2. 
Dit "O22 








IREF= 4 
JOH 6) SA oa 8. 
DS= 0.89 
DKG2= 2.70 











le - 2 
ITYPR= 2 
| CFORMI= 0.1514 
CFORM2Z= 0.1514 
| DS= 0.5 
DTI= 0.8 
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DECODE SUBROUTINE (PART 5) 










CC(1,1,1)=1.6E-03 
Car, 2,1) = 
CC(1,3,1)=0.0 


Os) 


Ce(2 71,71) =4:. 75-02 
CC(2,2,1)=-0.1097 


CC(2,3,1)=6.7E-02 


CC(3,1,1)=-0.1164 
CC(372,1)=-0.2176 


CC(3, 3,1) =-9. 7E-02 


CT(1,1)=0.834 
ER Bie Os er 


S <a (>) IBOTH-1 





{occ (1,1, IBOTH)= 3 





CGC, 3 ,1Bberi) =) 2 















COU, 6,1) = 






CC(2,2,1)= 






eet2 7,1) 
SE(S,2,1)= 
ee1 3, 3,1)= 
Che, 1 jel 
CT C23) =e 1 










Sal > 


6] 


CC(G@cl,l)= 
Cee). 2303 


RETURN 


> re 


(+ 


(0) 





7550 


CC (L727 r20TH)= -3.3478 
.8726 
Cees orH) —9tl. 2617 
GE(2,2,1BOTH)= -27.4162 
CC(2,3,IBOTH)= 21.7308 
| CC(3,1, IBOTH}= -23.7915 
CC(3,2,IBOTH)= 44.2607 
CC(3,3,IBOTH)= -14.4996 
CT(1,IBOTH)= 0.622 
CT(2,IBOTH)= 0.885 


- 1041 


- 1676 


CC (2,1,1)=-.1940 


aC MOM es: 


GCl27 2) =. 1646 
7.3E=-02 
ae Ob-O2 
2.44E-03 
aO32 

230 

























TRAJ SUBROUTIVE 
ENTRY 


\/ 


GE= ser ORM 
DM= DM1 
DKG= DKG1l 
MSTG= 1 
X= 0.0 

T= 020 
VX= VXA 
VY= VYA 
TH= FN 
YA= Y 

Y= ALT(IY) 


\/ 


(+) ; (=) 


/ 


q (0) 2 
D= DS+SL*U D= DMAX 


KUNGS 









* eats r = 


DTV= 1/G* (VY+SORT (VY**2+2.*G*Y) ) 
D= DTI 








I DNO 
ol gs 
.EQ.23 


(F) 





Seer TYPR= -1 IN FORTRAN 
TABLE]= O 


Sel PM 
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TRAJ SUBROUTINE (PART 2) 


ea 


DKG= DKG2 
DM2 


CFORM2 
FAS) 








DKG= DKG2 
DM= DM2 
CF= CFORM2Z 
TH= 0.0 


Oty 
“A DIV-D 


(0) 


ee 


RUNGE 


DTV=1/G* (VY+SORT (VY **2+2. *G*Y) ) 
T= T+DTV 

X= X+DTV*VX 

SWITCH= 0.0 













PLM TABLE1= 27 


4 
‘ MSTG= 6 


oO 


ip ; PLM TABLE 
MSTG= 11 MSTG= 0 





i 


‘ 








RUNGE SUBROUTINE 


ENTRY 


\/ 
GLOBAL 
VARIANCES 
van 


*} AD= A*D 


Yor Y 
VXO= VX 
VYO= VY 


/ 


RHO=2.4E-04-Y (6. 8E-08-Y (6. 7E-13) ) 





\/ 


Y=YO+AD* VY 
APl= AP2 
AN1= AN2 
VX@ VXO+AD*AN1 
VY= VYO+AD*AP1 









RHO=2. 4E-4-Y (6. 8E-8-Y (6. 7E-13) ) 


+e[T= THD | 


X= X+D* (VXO+AA* (VX-VXO) ) 
Y= YO+D* (VYO+AA*® (VY-VYO) } 
VX= VXO+D* (ANIL +AA*® (AN2-AN1) ) 
VY= VYO+D* (AP1L+AA* (AP2-AP1) ) 







il 


FA LS AY TN) PLM 





opr 1S TM IN PLM 


64 





Ot 
OA 
6 
OAA 


8 


9) 
OAA-AA aiid BO=S205 9 v 


€CO-ALE*ST .€T-EZL°9 
p LINW ans LTAWT 


onijans [T\ - 
‘ Ze, ) cox [* 
o-oo 
TdW 
CO-dLE°Z? 
S 
; - = [ere 


A 
; LTIIW 








INW TIOW 





CN 


XA 
OXA 


At Die. eternaneeli Med 2... 
(T Lewd) LUWHD SSHD0ud FONNY 





OA 


OA 
OAA-AA 
@) , Ln ILTIAW 
OHA 
—_ * Seemed A 
& aay LInW qav} ; * 
ee Tay 
AA 
ao" 
Gy a 
XA 
( a OXA 





rr 
(Z LuUWd) IYWHO SSHD0Ud GONNa 


OT 
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DERIV SUBROUTINE 


ENTRY 


GLOBAL 


VARIANCES 


V 


V= SORT (VX*VX+VY*VY) 
CM= V* (8.9E-044+3.2E~09*Y)+DM 


) (se) 


V 


IREG= 3 No IREG= 2 | 


vs 


oem OmGter (CC(IREG, 1 MSTG) +(CC(IREG, 2,MSTG) +CC(IREG, 3,MSTG) *CM) *CM) 


V 


HH= TH/V- RHO*CKDG*V 


AN2= HH*VX 
APZ— HHSVyY—G 


CP) So 


RETURN 





LTOW 


ans 


HH 


LTO | 


_————— 








LYOS Ga 
Tans 





INWHD SSdoodd ALYAC 
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APPENDIX B: 


Sa tibolL 
PL/M if 
Vaid ed 


FORTRAN 
A Al 
AA 

ALT 
AN1,AN2 
AP1,AP2 
GC 

CF 
CFORM1 
CFORM2 
CKDG 

CM 

er 

DEG 

DEL 

DKG 
DKG1 
DKG2 

DM 
DMAX 
DM1 

DM2 


PROGRAM DEFINITION OF VARIABLES 


DEF INTIIGN 
Runge Kutta parameter 
0.5/Runge Kutta parameter 
Weapon release altitude 
Runge Kutta variables (x direction) 
Runge Kutta variables (y direction) 
Viewer or dy; g0mcunvec COctLVeErents 
CKDG stretch factor 
GLWGastrerci sacuor Lor EnCe,Irst sage 
ChioUGustreteheractOnr FOL nme SecCOnds= tace 
Bomb coefficient times drag coerficient 
Mach number 
Moa ibcsOrmmNaewecUrve COcCii1c1ents 
Weapon release angle 
ArcTan (ejection velocity/weapon velocity) 
Shift in CKDG 
Simp eimemDG Lone tne fixst stage 
ShiteinwexDG Lor the, second stage 
Shift in Mach number 
FPaaeect In teeGacronmstep Size 
Shift in first stage Mach number 
Shift in second stage Mach number 


Dicceomats On estep Size 
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SYMBOL 
PiM if 


FORTRAN Varied 


DTI 


IDNO 
IREF 
TREG 
IV 
ry 
MSTG 


TH 
THETA 
U 

V 

VX 


DEFINITION 
Mit ciciesineeeration step size 
EqUdi~cOn tic =  imucStmrar tien! Inst 9s cape 


Parameter for Runte Kutta integration 
(value between zero and one) 


Acceleration due to gravity 
Total drag function 
Number of release angles integrated 


Key determining if two drag curves 
are required for the weapon 


Identification number of each weapon 
Reference for Mach curve coefficients 
Mach region index 

Number of velocities integrated 
Number of altitudes integrated 
Weapon stage index 

Conversion factor for radians 

Air density 

Ee 106. een ine variables 


SlppememOcnnrelent 1Or retarded snakeye 
PiiiconaccpLoyment time 


Elapsed time of flight from weapon release 
innwsty Masseot rocket 

Adjusted weapon release angle 

Mirerake: speed 

Velocity of weapon 


Velocity cemponent of V in the x 
dance tion 
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SYMBOL 
PL/M if 


FORTRAN Varied 


VXA 
VxXO 
VY 


VYA 
VYO 
VZ 


YO 
aT 


DEFINITION 


VX at weapon release 
VX at start of present integration step 


Velocity component of V in the y 
direction 


Naat Weapen release 

VY at start of present integration step 
Muzzle velocity of guns 

Weapon ground range from release 

Weapon aititude above sea level 

Y at start of present integration step 


Target alititude above sea level 





PORBRAN AND PL/M RESULTS 


IDNO VATS ALT DEG TOE | RANGE 


3000. a1. || i, 26!) wouls. 
3000. Si. |e lo, ai pscce 
3000. Eis || uO, S26 yaaa 
3000. “10. | 10.98) gies, 
3000. “1. | 8) 70es. 
3000. “100 )) 11, tel Gags 
3000. Si, | 10,30) Wales. 
3000. Sloe Oss 7292: 
3000. 10. 1h GG 2a) yee, 
3000. 10. ele 251) ys le 
3000. =102 |} 10.2515 fees) 
3000. =10, |) 1@ee4l Pease 
3000. Mew 22t 7547. 
3000. “NO |) a, Suh ee 
3000. 10, || Sieeten | yan 
3000. PHO 10626) 7516. 
3000. =i, | SOMES EVO 
3000. (iO. || 1625S) sh3- 
3000. -10. | 20.84} 4620. 
3000. 210 || TANS | Sale 
3000, 270. | WE, Sl) Bee 


3000. aoe ig oI oe 





3000. imu eG 7) 10040. 





FORTRAN AND PL/M RESULTS (continued) 
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